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METHOD AND SYSTEM FOR THE PREDICTION OF EARTHQUAKES 

The present invention relates to a method and a 
system for the prediction of earthquakes- In particu- 
5 lar, the present invention relates to a system and 
method for the detection and preventive localization 
of earthquakes; the system is installed in high seismic 
risk regions and consists of several receiving stations 
capable of detecting low frequency electromagnetic 
10 waves which are released from below the ground due to 
the mechanical stress of underground rock layers sub- 
jected to compression. 

The results of the detections are sent to a con- 
trol center capable of processing the data received to 
15 obtain the coordinates, an indication of the intensity 
and when the seismic shock in formation is expected 
(and the period of the following earthquake) . 

It is known that rocks which form surface layers 
of the earth's crust, when subjected to extremely 
20 strong mechanical compression stress, emit electromag- 
netic waves in the VLF band. 

The Applicant has observed that the measurement of 
these pre-seismic emissions received from special sta- 
tions distributed over a territory, and their subse- 
25 quent mathematic processing, allows the epicentre and 
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expected intensity of the seism to be identified with 
close approximation, and with the help of an accurate 
statistic study of previous shocks, it is possible to 
determine with sufficient precision the number of hours 
5 which preceding a destructive shock. 

The present invention is based on the fact that in 
seismic areas, the earth's crust is periodically sub- 
jected to strong tensions and compressions which, when 
reaching breakage, give rise to catastrophic seismic 

10 phenomena. 

It is known that seismic regions of the earth's 
crust mainly consist of rocks with high crystalline 
contents and that if these are subjected to strong me- 
chanical stress, they generate electricity (Piezoelec- 

15 tricity) , with the consequent generation of VLF elec- 
tromagnetic waves. These variations in electricity cre- 
ate electromagnetic waves with an intensity and fre- 
quency correlated to the mechanical stress which has 
generated them. 

20 These electromagnetic waves have frequencies rang- 

ing from a few tenths of Hertz to several tens of Kilo- 
herz; as the ground tends to attenuate more high fre- 
quencies, the deeper the phenomenon which has generated 
them, the lower will be the spectrum received. 

25 As it is known that the intensity of these elec- 
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tromagnetic emissions is proportional to the square of 
the forces applied per unit of rock volume, the Appli- 
cant has observed that, particularly after an accurate 
statistic study of the territory, it is possible to 
5 predict with a good approximation the intensity of the 
seism which is under formation. 

It should be noted that electromagnetic waves gen- 
erated in an area of mechanical anomaly, are attenuated 
with the square of the distance from this point to the 

10 receiving point, this means that the further away from 
the source, the intensity of the premonitory emissions 
of the seism will become increasingly weaker and con- 
cealed from disturbances generated by human beings, 
static discharges and foreign radio waves* 

15 For this reason, special analogical and digital 

(DSP) processing and filtering techniques of the sig- 
nals have been used in the present invention. 

The objective of the present invention is to pre- 
dict, on the basis of an analysis of premonitory VLF 

20 electromagnetic wave emissions, where the epicentre of 
the seism is located, the approximate intensity of the 
shock and when this will occur. 

In order to determine the epicentre of the seism, 
a detection network is used on the territory with re- 

25 ceiving stations situated at a distance of various tens 


-4- 

to approximately a hundred kilometers from each other 
(a minimum of three) , possibly in areas as far as pos- 
sible from sources of disturbance such as high voltage 
electric lines, cities, industrial centers. Each of 
5 these peripheral stations consists of a series of suit- 
able components housed in a metal container (protected 
or shelter) . 

These peripheral detection stations are programmed 
so as to be able to automatically send, via radio or 

10 telephone, the result of their detections to a central 
unit capable of processing the signals coming from said 
peripheral units, until the exact location of the epi- 
centre of the future earthquake is identified. 

An aspect of the present invention relates to a 

15 system for the prediction of earthquakes comprising at 
least three peripheral detection stations distributed 
in pre-determined positions over a territory, each sta- 
tion being capable of detecting the low frequency elec- 
tromagnetic waves which are released due to the me- 

20 chanical stress of the underground rock layers, a cen- 
tral processing station which receives signals relating 
to said electromagnetic waves from each of said periph- 
eral detection stations and processes said detections 
until an indication of the epicentre and entity of said 

25 earthquake is obtained, characterized in that each of 
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said peripheral stations comprises at least one ring 
aerial for each of the directions defined by a set of 
three Cartesian axes suitable for detecting said elec- 
tromagnetic waves, at least one of the ring aerials is 
5 an omnidirectional wide band aerial connected to a re- 
ceiver having an FFT process in real time suitable for 
detecting a preferred frequency bands for receiving 
said electromagnetic waves substantially without dis- 
turbance frequencies, this allows the reception of all 
10 ring aerials on said disturbance-free band free to be 
tuned in, thus widely increasing their sensitivity. 

A further aspect of the present invention relates 
to a method for the prediction of earthquakes, charac- 
terized in that it comprises the following phases: 

15 • detecting low frequency electromagnetic waves which 
are released due to the mechanical stress of under- 
ground rock layers, by at least three peripheral de- 
tection stations distributed in pre-determined posi- 
tions over a territory, 

20 • temporally and accurately marking the premonitory 
signals relating to said electromagnetic waves by 
means of a GPS receiver, 

• sending said marked signals to a central processing 
station, 

25 • calculating the differences in the receiving times 
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between a signal coming from a peripheral station and 
another signal coming from a different peripheral 
station. 

• calculating on the basis of said time differences of 
5 arrival of the signals, the coordinates of the epi- 
centre and depth of the earthquake by effecting a 
TDOA calculation. 

For the sake of greater clarity, the description 
of the invention refers to the enclosed drawings, which 
10 are provided for purely illustrative and non-limiting 
purposes, in which: 

- figure 1 schematically shows in its essential 
parts, the prediction system according to the present 
invention comprising at least three peripheral detec- 

15 tion stations and at least one central processing sta- 
tion; 

- figure 2 schematically shows one of the above 
peripheral detection stations of VLF electromagnetic 
waves according to the present invention; 

20 - figure 3 shows a block scheme of the operating 

components present in the above peripheral detection 
station according to the present invention; 

- figure 4 is a diagram which illustrates the pro- 
cedure with which a central monitoring unit processes 

25 the information received from the various peripheral 
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detection stations in order to locate the epicentre of 
the seism according to the invention; 

- figure 5 is a three-dimensional cross-section of 
a portion of the earth's crust where four peripheral 

5 detection stations belonging to the system according to 
the present invention, are operating. 

- figure 6 illustrates an amplifying circuit for 
ring aerials according to the present invention. 

With reference to the above figures, the predic- 

10 tion system according to the present invention com- 
prises at least three peripheral detection stations 2, 
2' and 2" distributed over the territory and at least 
one central processing station 3 which receives the de- 
tections effected by the peripheral stations. 

15 Each peripheral detection station comprises a sys- 

tem of three main balanced protected ring aerials 21, 
21' and 21" and at least one secondary wire aerial 
preferably underground. 

Of the protected ring aerials, there is preferably 

20 one for each direction from which the detection is ef- 
fected. For example, a first wide-band ring aerial 21 
arranged horizontally, is used, together with a second 
ring aerial 21' arranged vertically in a north-south 
direction, and a third ring aerial 21" arranged verti- 

25 cally, orthogonal to the previous one, arranged in an 
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east-west direction. In this way, detections can be ef- 
fected on the three directions forming a series of 
three Cartesian axes x, y and z. 

These aerials detect the electromagnetic waves 
5 which are generated from seismic phenomena as mentioned 
above . 

Weak signals leaving the magnetic ring aerials are 
amplified with special very low noise current amplifi- 
ers situated directly in the metallic boxes in the base 

10 of the aerials themselves in order to improve their 
disturbance rejection, the same container also contains 
the necessary circuits for the automatic tuning of the 
N-S E-W aerials, these signals are then sent to suit- 
able devices contained in a container 22, (after being 

15 adapted in impedance, amplified and filtered) where the 
signal is converted from analogical to digital. 

Said peripheral detection unit also comprises a 
GPS 23 receiver, at least one analogical digital con- 
verter for each aerial present in the unit 24, 24' and 

20 24", a local processing unit 25, a memorization device 
2 6, a telephonic modem 27 and a modem via earth or sat- 
ellite radio 28 equipped with an aerial 29 (for appli- 
cations in particularly inaccessible areas, a connec- 
tion via satellite can be used) . 

25 Once they are inside the container, the same sig- 
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nals detected by the above aerials are converted from 
analogical to digital, by means of said A/D devices. 

Once these signals have been converted to digital, 
they are sent to said central processing station after 
5 being sampled by means of a synchronization signal com- 
ing from said GPS 23 satellite receiver which also sup- 
plies the time-sample marks. 

Said peripheral detection stations are programmed 
so as to be able to automatically send, via earth or 

10 satellite radio modem or by said telephonic modem, the 
result of their detections to said central station ca- 
pable of processing, by means of a central processing 
unit, the signals coming from said peripheral units, 
until obtaining the exact localization of the epicentre 

15 of the future earthquake. 

The operating logic of the central processing cen- 
ter is schematized in figure 4, in which the processing 
of signals coming from four distinct peripheral detec- 
tion stations (Rx-A, Rx-B, Rx-C, Rx-D) is assumed. 

20 More specifically, this operation consists in 

processing the signals coming from the detection sta- 
tions, suitably equalized, and finding the correlation 
peak among them (the correlation consists in letting 
the signals received from the stations temporally pass 

25 numerically until, within a definite time lapse, they 
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coincide with each other) . 

The correlation peak serves to temporally align 
the signals received in the detection stations, thus 
being able to measure the relative delays [t(a-b), t (a- 
5 c) , t(a-d)] between the time marks of the signals re- 
ceived from the stations distributed over the terri- 
tory. 

In other words, the VLF signals received from the 
peripheral detection stations are marked with a time 

10 signal deriving from the GPS station, and subsequently 
the arrival time differences are calculated in the cen- 
tral processing station, by means of the correlation 
and GPS marking technique, as shown in figure 4, A pos- 
sible temporal inaccuracy of a marking with a GPS does 

15 not exceed a micro-second, thus providing a theoretical 
precision of .3 Km in the determination of the epicen- 
tre, assuming that this is. in a sole point. 

The delays are measured on several samples of sig- 
nals received from each detection station; the results 

2 0 are then averaged and, when a minimum datum variance is 
reached, the value is acquired. 

At this point, as the coordinates of the detection 
stations and receiving time differences of the premoni- 
tory signal are known, it is possible to obtain the co- 

25 ordinates of the epicentre and depth of the seism using 
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the known TDOA (Time Difference of Arrival) technique. 

This technique does in fact make it possible, with 
three detection stations, to operate in 2D i.e. the co- 
ordinates of the epicentre in x-y; whereas, with four 
5 or more detection stations it also ensures, in addition 
to a greater precision, the possibility of operating in 
3D, i.e. of also detecting the z axis which represents 
the depth of the epicentre. 

With respect to an assigned Cartesian reference 
10 system (x, y, z) as shown in figure 4, the incognito 
coordinates of the emission signal are obtained by 
solving the system: 

\Kx-x b Y + (y-y b y + (z-z b Y-il(x-x a y+(y-yJ + (z-*J=c 

• J(x-x c f +(y-yj +(z-0 2 - J(*-xJ +(y-y.f +(*-0 2 

Jix-xJ + {y-y d f + (z-z rf y -J(x-xJ +(y-yj +(z-zj ^ 

Wherein t(ab), t (ac) , t(ad) are the delays of the sig- 
nals received from stations B, C, D respectively with 
respect to the signal received in A; wherein the con- 
stant C is the propagation rate of the electromagnetic 
waves underground . 

With the Cartesian system indicated above, it is 
therefore possible to identify the coordinates (x, y) 
of the epicentre as well as its depth (z) . 


15 
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If there are only three peripheral detection sta- 
tions, the processing procedure should take into ac- 
count the delays of the signals of two stations with 
respect to the signals of the third; the relative cal- 
5 culation system would consequently consist of only two 
equations with two incognitos: specifically the coordi- 
nates (x, y) of the epicentre. 

According to the present invention, in order to 
improve the sensitivity of the overall receiving sys- 
10 tern, a method is used for better adapting and tuning 
said ring aerials as described below. 

It has already been mentioned that one of the ring 
aerials is a wide band aerial, i.e. an aerial capable 
of detecting frequencies for example of 10Hz up to tens 
15 of Khz. 

Each ring aerial is equipped with an aerial tuner 
capable of regulating the spectrum of frequencies which 
can be received from the aerials themselves. 

Said first horizontal ring aerial is a wide band 

20 aerial and is equipped with a microprocessor digital 
receiver (piloting electronic circuit) , which effects 
an FFT analysis of the signal received from the aerial. 
On the basis of the spectrum analyzed with said FFT, 
all the stable disturbance frequencies identified as 

25 disturbance are eliminated (in which the presence of 
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known disturbances is possible), (for example a and 
Omega signals used by submarines at 40 Hz and 80 Hz ap- 
proximately/ natural signals such as the discharge of 
streaks of lightning received on other bands and can- 
5 celled by the main signal, Shumann resonances at 7.7 
and 21 Hz, etc.). The remaining frequency bands which 
are substantially free of disturbances are selected for 
the reception of the electromagnetic waves of the pre- 
cursors. The piloting circuit analyzes and controls the 

10 aerial tuner in real time in order to allow a re-tuning 
of the aerials if a new disturbance is introduced into 
the system. The tuners of the remaining aerials are 
also preferably controlled by means of the FFT effected 
on the signal received from the wide band aerial. 

15 Another possibility for improving the sensitivity 

of the aerials is to provide them with a low noise am- 
plifying circuit 6 as illustrated in figure 6 which 
preferably uses an operational amplifier 61. 

The receiving ring aerial is accurately schema- 

20 tized and represented as an Ra resistance in series 
with an inductance La. Said aerial is equipped with a 
counter-reaction coil Lc having a counter-reaction re- 
sistance in series Rc. 

It is known that the sensitivity of an aerial con- 

25 nected to a current amplifier has two limits, with an 
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increase in the frequency , the impedance generated by 
the inductance La increases and consequently the sensi- 
tivity of the aerial decreases with higher frequencies, 
whereas the resistance Ra limits the general sensitiv- 
5 ity of the aerial. In order to avoid this, the reaction 
inductance Lc provides and additional La field, this is 
obtained by feeding this with the outgoing signal con- 
verted from voltage current (VCS) suitably regulated a 
regulation resistance R4, this opposes the inductance 
10 effect La bringing the total impedance of the aerial to 
almost zero, by extending its receiving frequency band 
upwards . 

The amplifying circuit comprises an RC circuit consist- 
ing of a CI condenser having a suitably predetermined 

15 capacity and a divider consisting of a resistance R2 
and an adjustable resistance R3, which when suitably 
regulated, provides a positive reaction which reduces 
the effect of Ra on all frequencies* 

The presence of a resistance R5 with a suitably 

20 calculated negative coefficient stabilizes the gain of 
the amplifier against the resistance variations of the 
copper wire of the aerial due to the variations in tem- 
perature . 

The amplifier is substantially a current amplifier 
25 in which a resistance Rl regulates the current gain, 
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and a condenser C2 stabilizes the operational 61 at ex- 
tremely high frequencies. The use of these techniques 
already known in electronics, applied in this particu- 
lar circuit, allows the efficiency and sensitivity of 
5 all the aerials used to be increased by numerous dB. 

Furthermore, in order to avoid all interferences, 
all the sensitive amplification and tuning circuits are 
arranged inside an airtight protected metallic con- 
tainer solidal with the aerial itself. 

10 According to the present invention, additional 

magnetic aerials can be used, capable of picking up ex- 
tremely low frequency signals, for example from 0.1 Hz 
to 10 Hz. These aerials comprise a cylindrical magnetic 
nucleus made of ferrite on which a copper wire is 

15 wound. The arrangement of the cylinder determines the 
direction x, y or z from which the aerial receives the 
electromagnetic signals . 

These aerials can also be advantageously equipped 
with said amplifying circuits. 

20 In order to improve the reliability of the system, 

each peripheral detection station comprises a radiogo- 
niometric localization device which is capable of de- 
tecting the direction of arrival in degrees of the 
electromagnetic signals. These directions are sent to 

25 the central processing station which effects a process- 
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ing, by crossing the directions of origin of the elec- 
tromagnetic waves in the various peripheral stations, 
to determine the epicentre of the seism. 

The point of origin of the signals is compared 
5 with the results obtained with the delay difference 
method and the radiogoniometric and mediate method, a 
map of the space distribution is then made in order to 
be able to track the borders of the generation area of 
the signals; this datum can give an indication on the 
10 intensity of the seism in formation; the vaster the 
area of origin of the signals, the stronger the seism 
will be- Studies effected indicate that an area of 20 
Km in diameter can give a seism of about the 7~ h Richter 
degree . 

15 A further estimate of the epicentre and entity of 

a seism can be obtained, combined with the electromag- 
netic wave system, by measuring the variations in the 
intensity and direction of the currents dispersed over 
the land by means of electric field measurers. In par- 

20 ticular, these extremely low frequency electric fields 
can be detected by means of NS-EW cross aerials con- 
sisting of metal wires arranged on the land, orthogonal 
to each other. The signals detected at the heads of 
these wires can be treated in substantially the same 

25 way as the signals coming from the electromagnetic 


waves , (marking with the GPS signal) and then subse- 
quently sent to the central processing station. 

Furthermore/ each peripheral station can substan- 
tially be in a rest position until one of the signals 
coming from any of the aerials exceeds a predetermined 
threshold. In reply to one of these events, the station 
can be activated and at the same time also allow the 
other peripheral stations to be activated through said 
central unit. 


